i>r 


SIMULATION  OF  SOLAR  KEATING  AND  COOLING  SYSTEMS 
USING  THE  CONTINUOUS  SYSTEM  MODELING  PROGRAM 


THO  CHING  HO 
.  S.,  Taiwan  University,  1973 


=l  MASTER'S  THESIS 


.tted  in  partial  fulfillme 


KANSAS  STATE  UNIVERSITY 


Approved  by: 


|<?7»  ACKNOWLKi)'>;:ENi':: 

C.  2      The   author  wishes    to  express   his    gratitude    to  his  major  advisors, 
Professors   R.    C.   Akins   and   B.    G.    Kyle,    for  their  guidance,    constant 
encouragement  and  valueable  advice    throughout   this    investigation  and 
preparation  of   this    thesis.      Appreciation   is   also  extended    to  Dr.    R.   L. 
Gorton  for  serving  on   the   author's    supervisory   committee  and    for   reviewing 
this    thesis   manuscript. 

for  their  continuous   encouragement   and    support  during    the   years   of 

Financial  support   provided  by   the  Engineering  Experimental  Station, 
Kansas   State  University,    is    gratefully  acknowledged. 


TABLE   OF   CONTENTS 

LIST  OF  FIGURES 

LIST  OF  TABLES 

NOMENCLATURE 

Chapter  1.   INTRODUCTION 

Chapter  2.   LITERATURE  REVIEW 

GENERAL  REVIEW 

THE  SIMULATION  METHOD 

CONCLUSIONS  OF  OTHERS  STUDIES 

Chapter  3.   SIMULATION  MODELS 

MODEL  A  -  SOLAR  ENERGY  HEATING  AND  CONVENTIONAL 
COOLING  SYSTEM 

THE  HOUSE 

ENERGY  COLLECTION 

ENERGY  STORAGE 

DOMESTIC  HOT  WATER  LOAD 

HOUSE  HEATING  LOAD 

HOUSE  COOLING  LOAD 

WEATHER  DATA 

AUXILIARY  ENERGY 

MODEL  B  -  SOLAR-ASSISTED  HEAT  PUMP  FOR  HEATING 
AND  COOLING 

ENERGY  STORAGE 

HEAT  PUMP 

AUXILIARY  ENERGY 

MODEL   C    -    SOLAR-ASSISTED    HEAT    PUMP   FOR  HEATING 
AND   HEAT   PUMP  CHILLED  WATER   FOR  COOLING 

ii 


CONTINUOUS  SYSTEM  MODELING  PROGRAM 
SIMULATION  RESULTS 
S  RESULTS 


Chapter  6.   COST  ANALYSIS 
Chapter  7.   CONCLUSION 
REFERENCES 
APPENDIX 

A.  SIMULATION  PROGRAM  OF  HEATING  FOR  MODEL  A 

B.  SIMULATION  PROGRAM  OF  HEATING  FOR  MODEL  B 

C.  SIMULATION  PROGRAM  OF  COOLING 

D.  WEATHER  DATA  DURING  THE  SIMULATION 

SOLAR  RADIATION   (OCT. -APR.) 
AMBIENT  TEMPERATURE   (OCT. -APR.) 
WIND  SPEED   (OCT. -APR. ) 


LIST  OF  FIGURES 


1  Solar  Energy  Heating  and  Conventional  Cooling  System. 

2  The  Flat  Plate  Sheet  and  Tube  Solar  Collector. 

3  The  Storage  Tank  Operation  For  House  Heating  and  For 
Hot  Water  Heating. 

4  Definition  of  Solar  Angles. 

5  Solar-Assisted  Heat  Pump  Heating  System. 

6  The  Heat  Pump  Heating  Cycle. 

7  The  Heat  Pump  Cooling  Cycle. 

8  Combined  the  Storage  Tank  and  the  Heat  Pump  for  Cooli 

1   Simulated  Cooling  Load  for  the  Hottest  Month. 

Lar  Energy  as 

3-5.10  Simulation  of  Auxiliary  Energy  Required  for  Hot 
Water  Heating. 

11-5.18  Simulation  of  Auxiliary  Energy  Required  for  Hous 


30  Values  of  CO. P.  for  MODEL  B  and  MODEL  C  during  the  Simulat 
1  Annual  Dollar  Savings  for  MODEL  A  and  MODEL  B. 


1  Design  Parameters  of  the  Solar  Collec 

2  Design  Parameters  of  the  Storage  Tank 


of  Solar  Declinat 


and  MODEL  B. 

5  Comparison  Between  the  Simulated  Monthly  Performance 
and  the  Estimated  Monthly  Performance. 

6  Values  of  the  Heat  Pump  CO.  P.  for  MODEL  B  and  MODEL  C 
in  the  Averaged  Summer  Weather  Condition  in  Manhattan. 

1  Annual  Dollar  Savings  for  MODEL  A  and  MODEL  B„for 
Collector  Areas  of  400,  600,  800,  and  1000  ft  at 
Power  Rates  of  $0.02,  0.03,  0.04,  and  0.05/KWh. 


NOMENCLATURE 


L   outside   surface,    ft 


■a,    ft 
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.  of  wall,  ft2 
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bond  conductance,  Btu  hr^ft"10?"1   (see  Figure  3.2) 
.e  capital  cost  per  unit  area  of  collector,  $/ft2 
le  annual  cost  of  the  conventional  system,  $ 

le  capital  cost  of  additional  equipment  for  a  solar  system,  $ 
.e  capital  cost  of  an  electric  furnace,  S 

le  capital  cost  of  a  heat  pump,  $ 

le  annual  cost  of  maintenance,  $ 

le  annual  cost  for  MODEL  A,  $ 

ie  annual  cost  for  MODEL  B,  $ 

.efficient  of  performance,  defined  in  Eq.  (47) 

leal  coefficient  of  performance,  defined  in  Eq.(46) 

ie  unit  cost  of  power,  $/KWh 

ie  capital  cost  of  storage (medium,  container,  and  insulation),  S 


tube  diameter,  ft 

the  annual  dollar  savings  for  MODEL  A,  S 

the  annual  dollar  savings  for  MODEL  B,  $ 

a  function  of  h  and  N,  defined  in  Eq.(lO) 

fin  efficiency,  defined  in  Eq.(4) 

collector  efficiency  factor,  defined  in  Eq.(3) 

>r  heat  removal  factor,  defined  in  Eq.(2) 

insionless  ratio  of  the  solar  heat  gains  to  the  incident 


ficient  between  the  fluid  and  t 


led    in  Eq. 

ri4i 

ied    in  Eq. 

(15) 

led    in  Eq. 

fib) 

led    in  Eq. 

(61) 

led    in  Eq. 

(62) 

transfer 

coef 

-ft-oF- 


.-1 


J  heat  transfer  coefficient,  Btu  hr   ft  *"UF 

:  gain  due  to  infiltration  and  ventilation,  Btu  hr" 
:  gain  through  the  roof,  Btu  hr" 

:  gain  through  a  wall  due  to  air  temperature  diffei 


=at  gain  through  a  wall,  Btu  hr 
:tual  intensity  of  the  solar  radi; 


actual  solar  radiation  impinging  on  the  collector,  Btu  hr 

measured  solar  radiation  on  horizontal  plain,  Btu  hr"Xft" 

annual  factor  of  investment 

fin  thermal  conductivity,  Btu  hr^ft'V1 

local  latitude,  39°  for  Manhattan 

defined  in  Eq.(5) 

collector  fluid  capacitance  ra 

cooling  load  fluid  capacitance 

house  heating  fluid  capacitanc 

capacity  of  the  storage  tank  section  1,  Btu  °F 


--v-i 


nr-V*-1 


of  t 


i  3, 


of  the  absoi 


,-1 


Met 


power  requirements  for  heating  for  MODEL  A,  KWh 

power  requirements  for  cooling  for  MODEL  B,  KWh 

power  requirements  for  heating  for  MODEL  B,  KWh 

power  requirements  for  cooling  for  the  conventional 


.  for  h 


5  for  t 


lourly  auxiliary  energy  required  for  hot  water  heating,  Btu  h 
Long-term  cumulative  auxiliary  energy  required  for  hot  water 


hourly  auxiliary  energy  required  for  house  heating,  Btu  h 

hot  water  heating  load  supplied  by  the  solar  system,  Btu 
the  house  heating  load,  Btu  hr" 


ambient  temperature,  °R 
temperature  of  the  basement,  °F 

operation,  °F 


;ign  temperature,  UF 

:e  of  the  fluid  at  the  collector  inlet,  °F 

re  of  the  fluid  at  the  collector  outlet,  °F 


temperature  of  i 


sturning  from  the  heating  load, 


■  in  section  1  of  the  storage  tank,  °F 
in  section  2  of  the  storage  tank,  °F 


heat  pump  condensing  temperature,   R 

house  equivalent  heat  loss  coefficient,  Btu  hr_1ft" 

heat  loss  coefficient  of  the  storage  tank  section  1 


aefficient  of  the  s 


heat  loss  coefficient  of  the  roof,  Btu  hr"1ft"2°F"1 
heat  loss  coefficient  of  a  window,  Btu  hr"1ft"2°F"1 
heat  loss  coefficient  of  a  wall,  Btu  hr"1ft"2°F"1 
wind  s  eed  mil  h  "1 

the  distance  between  the  tubes,  ft 

solar  attitude  above  the  horizontal,  degrees 

solar  declination,  degrees 

fin  thickness,  ft 

collector  tilted  angle  from  horizontal,  degrees 

emittance  of  plate 


le   incident   angle,    degrees 

le    incident   angle   on   the  49°   tilted   collector   surfac 
tefan-Boltzmann  constant,    0.1712   x  10"      Btu  hr"1ft"'! 


CHAPTER  1.   INTRODUCTION 

have  shown   us    that    the  world's    fossil    fuel   resources,   which  might  be   con- 
sidered a    form  of    stored    solar   energy,    will    soon  be  depleted    in  meeting  man' 

the  supply   of   oil  will   fail    to  meet    increasing  demand  before    the  year    2000, 
most   probably  between   1985   and   1995,    even   if   energy   prices   rise   50%  above 
current   levels    in  real    terms.      The    task   for    the  world  will  be   to  manage  a 
transition   from  dependence   on   oil    to   greater    reliance  on  other   fossil 

Evidently,    the   renewable  resources,    especially   solar,    are   likely    to 

ceptions   lead    us   naturally    to   consider    the  ways    of   collecting,    storing, 
and    using  solar    energy. 

The  use  of   solar    energy    is   not   a   new   revelation.      Since   time   immemorial 
priests    ignited    temple   fires  by   concentrating   solar   energy,    thus   adding 

in   solar    energy   history  was    in   1772  when   a  French   chemist,    Lavoisier, 
generated    temperatures   up   to   1, 750°F  by   using    two   hollow  glass   lenses 
filled  with  white  wine   to   focus   sunlight.      Then   in  1870's,    John  Ericsson 
developed   solar   engines.      He   found    that    they    could   not   compete  with   fuels 
of    that    time.      After    that,    in    the   late   19th   century   and   early    20th    century, 
many    solar-assisted    instruments  were   designed   and   constructed.      These 


In  recent  years,  research  and  development  on  direct  use  of  solar  energy 
has  led  to  a  much  better  understanding  of  the  principles  of  construction 
of  solar  systems.   These  solar  systems  might  be  classified  as:   (A)  Solar 
Energy  for  Water  Heating  and  Steam  Generation;  (B)  Solar  Energy  for 
Desalination  and  Drying;  (C)  Solar  Energy  for  Building  Heating  and  Cooling; 
and  (D)  Direct  Solar  Production  of  Electricity. 

and  hot  water  is  the  most  advanced  and  the  most  competitive  with  fuels. 
Practical  systems  are  available  which  require  the  use  of  glass-covered 
panels  on  part  of  the  roof  and  a  small  space  for  overnight  heat  storage, 
usually  in  a  basement.   An  auxiliary  furnace  supplies  heat  when  necessary  - 

There  are  now  upwards  of  200  solar  heated  buildings  in  the  United  States. 
A  few  of  them  are  well-engineered  and  instrumented  systems,  primarily  in 
universities  and  research  institutes,  such  as  the  Colorado  State  University 
Solar  House  I,  II,  and  III,  the  MIT  House  IV,  and  etc..   Comparison  with  the 
present  costs  of  conventional  energy  shows  that,  solar  heating  is  not  competi- 
tive with  natural  gas  at  today's  prices,  but  it  is  considerably  cheaper  than 

development  are  still  needed  to  make  them  better.   In  this  thesis,  a  research 
result  is  reported  on  Solar  Energy  for  Building  Heating  and  Cooling.   Three 
models  of  solar  heating  and  cooling  systems  were  simulated  on  a  digital  com- 
puter using  C.S.M.P.  (Continuous  System  Modeling  Program).   MODEL  A,  called 


SOLAR  ENERGY  HEATING  AND  CONVENTIONAL  COOLING  SYSTEM,  is  mainly  composed  of 

electric  furnace  to  supplement  solar  energy  in  the  winter,  a  conventional 
air  conditioner  for  air-conditioning  in  the  summer,  and  the  necessary  pumps, 
piping,  and  controls.   MODEL  B,  called  SOLAR-ASSISTED  HEAT  PUMP  FOR  HEATING 
AND  COOLING,  is  similar  to  MODEL  A,  with  the  exception  that  a  heat  pump  is 

also  used  for  house  heating.  MODEL  C,  called  SOLAR-ASSISTED  HEAT  PUMP  FOR 
HEATING  AND  HEAT  PUMP  CHILLED  WATER  FOR  COOLING,  is  equipped  all  the  same 
as  MODEL  B,  but  the  ways  of  operating  the  heat  pump  for  house  cooling  are 
different. 

Manhattan,  Kansas.   The  hour-by-hour  performance  of  each  model  was  studied 
using  the  19  76-19  77  hourly  weather  bureau  data  for  the  Manhattan  area. 
Several  collector  areas  were  tested  in  these  models  and  a  cost  analysis 
was  used  to  determine  the  minimum  cost  system.   The  results  also  indicate, 
under  what  conditions  in  Manhattan,  which  of  the  solar  energy  systems  could 
be  economically  competitive  with  conventional  systems. 


CHAPTER   2.      LITERATURE  REVIEW 

GENERAL   REVIEW 

Space  heating  with  solar   energy,    after  about   35   years   of   limited 
development   effort,    has   recently  become  a   commercial   reality,    and    it    is 
likely    to   grow  into  a  major   energy   supply  and   industrial  activity.      The 
most   common   type   of   system  for   space  heating  with   solar   energy   involves 
heat   collection   in   a   liquid   passing    through   a   flat-plate   collector   and 

The   trend   in   liquid    collector   design  and   construction   is    toward  a 
factory-built  metal  box   containing   a  metal   absorber   plate   coated  with 
black  paint  or   a   selective  radiation   surface,    and   provided  with   insulation 
beneath    the  metal  absorber.      Tubes    in    the  absorber   plate  provide  passages 

restrict  heat   escape.      The  assembly  of   individual   collector  modules    into 

between   adjacent   collector   panels. 

Collectors   require   several    types   of   materials,    the  absorber  plate  being 

thin   oxide   film  on   a  bright  metal   substrate,    and    the   glazing  being   a  durable, 
heat  resistant  plastic  or   glass,    usually    tempered. 

Throughout  most  of    the  United    States,    freezing  may  occur   on  winter 
nights,    so  water  must  not   remain   in   a    typical   collector  when  no   heating 
is    taking   place.      One  method    for   dealing  with    this   problem   is   an  automatic 
draining   of    the    collector  whenever    the   circulating   pump    is   not   operating. 
The  other  method    is   by   use  of  a  non-freezing   liquid    in    the   collector. 

Transfer   of  heat   from   the   storage    tank   to    the  building   is   accomplished 
either  by   circulation   of   heated  water    through   radiators   or   convectors    in 


the   rooms,    or  by   its   circulation    through  a    fan   coil    exchanger   in   the  duct 

Although  not  as  widely  utilized    as   liquid    types,    a   few   solar   air  heating 

simplest   form,    a   solar   air   collector   closely  resembles    the  liquid   heating 
type.      The   essential  difference   is    the   circulation  of   air   in  contact  with, 
and   usually  beneath,    the  black  absorber   plate.      Rather    than  being    in    tubes, 
the   air  may   contact    the   entire   surface  of    the   absorber   plate.      To   enhance 
the   heat    transfer   coefficient,    fins,    corrugations,    or  other    types   of    extended 
surfaces   on    the  absorber  may  be   employed. 

Heat   storage   in  an  air   system,    for    the  best   collection  efficiency,    should 
also   extract   substantially   all    the  useful  heat   from   the  hot   air   stream   so 

enced  with   air   returning   from  the   rooms.      A  highly   stratified    storage 
capability   is   therefore  desirable.      The    thermal   properties   of   loose   solids, 

The  very   large   surface  area   of    the   pebbles    (commonly   2    to    3   cm  in  diameter) 
and    the  high  porosity   of   a  bed   of   uniformly   sized   solids   results    in   rapid 
transfer   of   heat. 

Supply   of   heat    to    the  rooms    from  the   storage  unit    is   easily   accomplished 
by  circulating   room  air    through   the  pebble-bed    in   a   direction  opposite    to 
that   employed   in    the   storing   cycle.      Room  air    is    thus  heated  by   contact 

Although   long-term   performance  data   on   liquid    systems    for   space  heating 
are  not  available,    certain   comparisons   can  be  drawn  between   the  liquid   and 


of  steel  or  aluminum  occurs.   There  is  no  possibility  of  freezing  the 
collector  fluid,  so  decomposable  antifreeze  compounds  or  other  organic 
liquids  are  not  required.   Boiling  can  not  occur,  so  replacement  of  collector 
fluid  is  never  needed.   Fluid  leakage,  although  not  precluded,  involves 
minimal  corrective  cost  and  no  secondary  damage.   Some  long-term  operating 
experience  with  an  air  system  shows  that  maintenance  is  not  a  significant 

to  that  of  the  building  itself. 

Although  the  basic  materials  in  the  air  system  may  be  cheaper  than 
those  in  the  liquid  type,  more  space  in  the  building  is  required  for  the 
heat  storage  units  and  connecting  ducts.   Interconnecting  conduits  and 

Lof,  et  al.  [4]  compared  the  performance  of  solar  heating  with  air  and 
liquid  systems  and  found  that  under  equal  conditions,  both  the  liquid  and 
the  air  solar  heating  systems  can  provide  space  heating  (and  expected 
service  water  heating)  with  nearly  equal  performance.   They  concluded  that 
an  air  system  should  have  more  appeal  for  residential  space  heating  where 

solar  heating  system  may  be  favored  in  larger  commercial  space  heating 
installations  where  maintenance  is  available  and  space  for  heat  storage 
and  equipment  conduits  is  of  high  value. 

The  addition  of  solar  cooling  to  heating  can  evidently  improve  the 

systems;  and  solar-related  systems  which  are  not  solar  operated,  but  which 


Operation  of  absorption  air  conditioners  with  energy  from  flat-plate 
collector  and  storage  systems  is  the  most  common  approach  to  solar  cooling 
today.   Heat-operated  cooling  cycles,  based  on  absorption  of  refrigerant 
in  liquid  absorbent  solutions  or  absorption  of  refrigerant  on  solid  adsor- 

step.   These  may  be  closed  cycle,  e.g.  the  LiBr-H20  or  NH3-H,0  cycles,  or 

atmosphere. 

A  solar-mechanical  cooling  system  that  has  received  some  attention  in 
recent  years  couples  a  solar-powered  Rankine  cylcle  engine  with  a  more 

with  this  system  are  basically  the  problems  associated  with  generating 
mechanical  energy  from  solar  energy.  The  storage  tank  transfers  energy 
to  a  heat  exchanger  which,  in  turn,  transfers  energy  to  a  heat  engine.   The 

When  solar  energy  equipment  is  installed  for  the  purpose  of  heating  a 
building,  some  of  this  equipment  can  be  used  to  cool  the  building,  but 
without  the  direct  use  of  solar  energy.   These  systems  are  classified  as 
solar-related  air  conditioning  systems,  such  as  heat  pump  systems,  sky 

Combining  the  solar  system  with  the  heat  pump  for  heating  and  cooling 

the  heat  pump  can  extract  heat  from  the  storage  tank  at  a  lower  water 
resulting  in  a  higher  collector  efficiencies. 


In  a  solar-heat  pump  system  which  supplies  heating  for  a  residence, 
the  heat  pump  can  be  used  to  supply  cooling  during  the  summer.   In  this 
case,  cooling  is  supplied  in  a  conventional  manner  and  solar  energy  is  not 
utilized.   (During  the  summer  solar  energy  in  such  a  system  will  probably 

systems  can  be  considered  in  three  distinct  configurations.  In  the  in-line 
mode  the  heat  pump  is  between  the  storage  tank  and  the  load  and  all  energy 
supplied  to  the  heat  pump  comes  from  the  solar  storage  tank.   In  the  parallel 

the  heat  pump,  which  uses  ambient  air  as  the  source,  acts  as  the  auxiliary 
for  the  solar  energy  system.  In  the  dual  source  mode,  the  source  of  energy 
for  the  heat  pump  can  be  either  the  solar  heated  tank  or  ambient  air  depending 

It  is  clear  that  there  is  a  wide  spectrum  of  potentially  interesting 
methods  for  using  solar  energy  for  building  heating  and  air-conditioning, 

extensively  studied  while  some  are  still  completely  undeveloped.   Many 

systems,  such  as  TRNSYS,  SOLAR,  NECAP,  MINI-SHAC,  CERES,  and  CAP  [6].   The 

while  the  latter  three  were  developed  to  supplement  the  first  three  by 

These  programs  make  the  analysis  of  solar  heating  and  cooling  system  easier 
and  quicker. 

The  future  of  solar  heating  appears  assured  by  the  growing  scarcity  of 
fossil  fuels  and  by  continued  improvements  in  solar  technology.   The  future 


one  alone. 

THE  SIMULATION  METHOD 

There  have  been  two  approaches  to  the  design  of  solar  energy  systems. 

term  system  performance  from  suitable  averaged  meteorological  data.   The 
second  and  more  popular  approach  has  been  to  formulate  a  computer  simulation 

The  transient  behavior  of  solar  energy  systems  using  mathematical  simula- 
tion models  has  been  studied  since  1967.   The  technique  is  to  write  an 
energy  balance  equation  for  each  component  of  the  system.   In  general,  the 

the  temperature  of  the  fluid  entering  the  collector.   This  in  turn,  for 
many  systems,  is  the  same  as  (or  is  closely  related  to)  the  temperature  in 
the  exit  portion  of  the  storage  unit.   The  outlet  temperature  from  the 
collector  becomes  the  inlet  temperature  to  the  storage  unit.   In  these 
equations,  time  is  the  independent  variable  and  the  solution  is  in  the  form 
of  temperature  as  a  function  of  time.   In  addition,  it  is  possible  to  inte- 
grate energy  quantities  over  time  and  thus  assess  the  energy  distributions 
for  the  system. 

tion,  the  amount  of  energy  delivered  from  the  solar  collector  to  meet  a  load 
and  the  quantity  of  auxiliary  energy  required.   The  simulation  also  indicates 


whether  the  temperature  variations  for  a  particular  system  design  are 
reasonable,  for  example,  whether  a  collector  temperature  would  rise  above 
the  boiling  point  of  the  liquid  being  heated. 

The  ability  to  predict  dynamic  behavior  and  make  energy  balances  over 

mating  the  effects  of  design-variable  changes  on  system  performance;  these 
design  variables  might  for  example,  include  selectivity  of  the  absorbing 
surface,  number  of  covers  on  the  collector,  collector  area,  and  so  on. 

these  systems  function  and  can  be  a  guide  for  experimentation. 

Computational  simplicity  in  a  simulation  model  is  needed  to  allow  an 

of  climates  at  reasonable  computing  cost.   Minor  details  in  the  component 
models,  from  which  the  system  model  is  formulated,  have  only  a  small  effect 
upon  the  long  term  system  performance  and  should  be  excluded.   A  properly 
formulated  simulation  model  provides  an  estimation  of  the  performance  of  a 

CONCLUSIONS  OF  OTHERS  STUDIES 

In  this  section,  some  major  conclusions  which  have  resulted  from  other 
studies  of  solar  energy  systems  are  reviewed.   Sheridan,  Bullock,  and  Duffie 
[7]  were  among  the  first  to  study  the  transient  behavior  of  solar  energy 
systems  using  mathematical  simulation  methods.   The  only  major  conclusion  of 

heating,  Close  [8]  presented  a  method  which  he  contends  is  useful  in  the 
analysis  and  design  of  all  transient  systems  of  this  type.   His  approach 


begins  with  the  formation  of  a  mathematical  model.   The  model  is  cast  i 

identified.  Then  factorial  analysis,  a  statistical  design  technique,  i 
used  to  determine  the  effect  of  these  design  variables  upon  the  system 
performance  or  economics  from  simulation  studies  with  a  digital  compute 


He  sbowe 


achieved  by  dividing 


fully  mixed.  He  compared  the  results  for  a  three-section  tank  and  a  six- 
section  tank  and  concluded  that  only  small  differences  exist  for  the  two 
models. 

Tybout   and   Lof    [9,10]    developed   a   digital   computer   simulation  program 
of  a   solar  heating  system;    they   used   a  model   of   collector   performance  and 
an   estimation   of    the   space  and  water   heating  requirements    to   calculate  hourly 
energy  balances    for   a   fully-mixed    fluid    storage   tank.      The   space  heating 

cost   of   a    solar   heating   system  was    estimated   by   calculating   the  yearly   cost 

the   estimated    initial    investment. 

Despite   the   simplicity  of    their  model,    the   economic   study  of  LSf   and 
Tybout    [10]    has   resulted    in  a  number   of    important   conclusions,   most   of  which 
are   in  agreement  with   the  conclusions   of   other.      Their   conclusions   are: 

(1)      The   effects    the  collector    thermal   capacitance  upon   the  overall 
system  performance  are  negligible  for   all    common   flat-plate 

plate  surfaces  were  optimum   in   even   the   coldest   climate. 


(3)   The  tilt  of  the  collector  from  horizontal  for  most  economical 

operation  of  solar  heating  systems  in  winter  should  be  the  latitude 
plus  ten  to  twenty  degrees. 

heating  to  improve  the  economics  as  well  as  to  ensure  reliability. 

system. 
(5)   The  optimum  energy  storage  capacity  of  solar  heating  systems  is 
the  thermal  equivalent  of  ten  to  fifteen  pounds  of  stored  water 
per  square  foot  of  collector  area. 
Based  on  1970  fuel  costs  and  very  optimistic  collector  costs  (S2/f t2) , 
they  [10]  also  concluded  that  solar  heating  may  be  competitive  with  conven- 
tional heating  systems;  their  simulations  indicated  that  the  cost  of  solar 

but  less  than  with  electricity  in  six  of  the  eight  locations  examined. 

The  long-term  performance  of  the  Colorado  Solar  House  has  been  reported 

by  Lof  [11],  has  been  operating  continuously  since  1957  with  no  maintenance. 
The  house  has  296  m2  of  living  area.   The  collector  is  an  Overlapped-Glass 
Plate  Solar-Air  Heating  having  a  cover  glass  area  of  49.1  m2 .   Solar  heat 
is  stored  in  a  rock  bed  of  primarily  granitic  rock  approximately  1.3  -  2.5 
cm  in  diameter.   The  results  showed  that  approximately  20%  of  the  total 
heat  load  was  from  solar  energy. 

Klein,  Beckman,  and  Duffie  [12]  designed  a  solar  space  and  water  heating 
system  for  residences.   The  system  considered  uses  liquids,  either  water  or 
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collector  is  used  to  transform  incident  solar  radiation  into  thermal  energy. 
This  energy  is  stored  in  the  form  of  sensible  energy  and  used  as  needed  to 
supply  the  space  and  water  heating  load.   In  their  simulation  study,  the 

of  Liu  and  Jordan  [14].   A  simple  space  heating  load  model,  the  energy  per 
degree-time  model  (to  be  discussed  in  Chapter  3),  was  used  to  estimate  the 
heating  demand.   The  simulation  was  using  TRNSYS  [15],  a  general  simulation 
program  for  solar  energy  systems.   Based  on  the  collector  cost  of  $80/m  , 
they  concluded  that  the  fuel  costs  must  be  about  $6/106KJ  ($6 .33/106Btu) 
before  a  dollar  savings  using  solar  heating  can  occur. 

Hunn,  et  al.  [16]  developed  a  model  which  simulates  the  effects  of 
hourly  weather  conditions  on  the  performance  and  cost  of  a  combined  solar/ 
conventional  heating  system  for  building  in  cold,  cloudy  climates.   The 

climates  is  economically  competitive  with  fuel  oil  heating  only  if  the  price 
of  oil  doubles  to  approximately  80c/gal. 

Lof  and  Tybout  [17].   The  mathematical  model  of  their  systems  was  modified 
to  include  heat-operated  absorption  cooling  units  of  the  aqueous  lithium 
bromide  type.   Approximately  one  hundred  analyses  of  hourly  cooling  and 
heating  performance  for  a  full  year  in  eight  locations  were  made.   They 

(including  water  heating)  offers  lower  costs  than  does  either  application 
alone.   Optimum  storage  for  heating  and  for  combined  use  is  about  50  Kg 
per  m2  of  collector  area. 
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Another  simulation  study  of  solar  heating  and  cooling  system  was 
accomplished  by  Butz,  et  al.  [18].   Their  system  includes  a  water  heating 

source,  and  associated  controls.  The  simulation  was  made  in  the  climate 
of  Albuquerque,  New  Mexico.  The  results  showed  that  this  system  would  be 
competitive  with  conventional  means  of  heating  and  cooling. 

The  early  analysis  by  Jordan  and  Threlkeld  [19,20,21]  and  Lof  [22]  showed 
that  solar  heat  pump  systems  are  economically  feasible  through  most  of  the 
United  States.   The  General  Electric  Phase  0  report  [23]  extended  the  study 

locations  with  a  variety  of  different  heating  and  cooling  modes.  A  signi- 
ficant advantage  was  shown  for  systems  with  dual  evaporators  that  switched 

Simulations  of  the  detailed  dynamic  behavior  of  solar-augmented  heat 
pump  systems  have  been  made  by  Freeman,  et  al .  [24]  and  Bosio  and  Suryanarayan 
[25].   From  these  and  other  studies  [26,27],  a  few  basic  solar  heat  pump 

advantageous.   These  systems  have  the  following  features  in  common:   an 

evaporator  side,  the  capability  to  utilize  ambient  air  as  the  heat  pump 
source  when  storage  is  depleted,  the  capability  to  bypass  the  heat  pump 
when  storage  temperature  is  high  enough,  rejection  of  energy  to  ambient 
air  in  the  cooling  mode,  "second-stage"  auxiliary  space  heating  by  duct 

for  domestic  hot  water  using  "unassisted"  solar  energy  -  i.e.,  collected 
solar  energy  which  is  not  augumented  by  the  heat  pump. 
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Many  other  simulation  results  for  solar  heating  and  cooling  systems 
have  been  reported  in  the  literature.  Most  of  then  have  similar  conclusions 
to  those  mentioned  above  so  have  not  been  reviewed.  Based  on  some  of  these 
important  conclusions,  three  solar  heating  and  cooling  systems  were  modeled 
and  simulated  in  this  study.   They  are  discussed  in  the  next  chapter. 


CHAPTER    3.      SIMULATION  MODELS 
A   feasibility   study   of    solar   energy  heating  and   cooling  systems  was 

Kansas.      Three  solar  space  heating  and   cooling  systems    (including  hot  water 
service)    which   use  water  as    the   energy    collecting  and   storing  medium,  were 
modeled,    examined,    and  compared   in    this   study.      In  order    to  make   the  pro- 
gramming  easier  and   quicker,    the   simulation  programs  were  written   using 
C.S.M.P.    (Continuous   System  Modeling  Program).      In  writing   the  programs, 
the  connecting  pipes  were  assumed  well-insulated,    with    no  heat    loss    through 

A  typical  two  story  house  located  at  Manhattan,  Kansas  was  chosen  as 
the  simulations  basis.  The  hourly  weather  bureau  data  of  1976-1977  were 
used  in  the  simulations.  The  three  models  were:  MODEL  A  -  SOLAR  ENERGY 
HEATING  AND  CONVENTIONAL  COOLING  SYSTEM;  MODEL  B  -  SOLAR-ASSISTED  HEAT 
PUMP  FOR  HEATING  AND  COOLING;  and  MODEL  C  -  SOLAR-ASSISTED  HEAT  PUMP  FOR 
HEATING  AND  HEAT  PUMP   CHILLED  WATER  FOR   COOLING. 

MODEL  A  -    SOLAR   ENERGY   HEATING  AND   CONVENTIONAL   COOLING   SYSTEM 

A   schematic  diagram  of    the  MODEL  A   is    shown   in  Figure   3.1.      A   solar 

were   equipped    for  house  and  hot  water  heating.      A  conventional   air   condi- 
concerned  with    the  weather,    the   energy   collection,    the  collected    energy 
in   the   following  pages. 


Figure  3.1     Solar  Energy  H 


THE   HOUSE 

A  well-insulated,    40  x   25  x  18   ft    (length  x  width   x  height),    two   story 
house  with  main   door    facing   north  was   selected    for   the   simulations.      The 

were  assumed    to  be   20%  of    the   total   side  wall   area.      The  basement  was 

ENERGY   COLLECTION 

The   flat-plate   collector,    the   collector   pump,    and    the  control   system 
are   involved    in   collecting   solar   energy  and    transferring   energy    to    the 

as    shown  in  Figure   3.2,    are:      the   "black"   solar   energy   absorbing   surface, 
with  means    for    transferring    the  absorbed   energy   to  a   fluid;    envelopes    trans- 
parent   to    solar   radiation   over    the   solar  absorber   surface  which   reduce 
convection  and   radiation   losses    to    the  atmosphere;    and   back  insulation    to 

The   solar   collector  was  modeled  based   on    the  information   from   the  book 
of  Duffie  and  Beckman    [28].      As   recommended  by    them,    the   thermal   capacitancs 
of   the   collector  was   neglected.      The  rate  of   energy   collection,    Q    ,    of  a 
flat-plate  collector  was    expressed  as: 

Qu  =  AcFR[Ic(TC°    "  VTf,l  "   Ta)]    =    (lil  Vc    (Tf,o  ~   Tf    i}  (1) 

The  collector  heat   removal    factor   F   ,   which   relates    the  actual   useful   energj 
gain  of  a   collector    to    the  useful    energy   gain   if    the  whole   collector   surface 
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Lgure  3.2  The  Flat  Plate  Sheet  and  Tube  Solar  Collec 


Where  F',    the   collector   efficiency   factor,    which   represents    the  ratio  of 
the  useful    energy   gain    to    the  useful    energy    gain   if    the   collector  absorbir 


F'    =  - 


_1_ 


lUL[W-(W-D)F]         C  TiD.hf    .    ' 

The  fin   efficiency,    F,   was    defined  as: 

In   the  above   equation,    Eq .    (4),    m   is   defined  as: 

- "  &**  ») 

In  order    to  apply   Eq .    (1)    to    the   simulation  work,    the   collector  overall 
energy   loss   coefficient,    U"L,    needs    to  be   calculated   first.      The   following 
equation   for   L'L  was  as   developed  by   Klein    [29]: 

,  o(T'2  +  T'2)    (T'+T') 


0,    - 


L        3.829  N 1_  +  1 -    2N+f-l 


T'   -   T' 


1.33  h  e  +0.05N(1-e 


P  g 


i  for  calculating  Tn>  hw>  C,  and  f  were  as  follows 
Q  /A 

ti[ 

hw  -   1.004  +  0.3  V 

C     =   365.9    (1  -    0.00883E  +  0. 0001298s2) 

f      =    (1  -   0.227h    +  0.016h    2)     (1  +  0.091N) 


Lmulation  equations  vv 


VI:,:- 


r  energy  by  cover 


affects  losses  from  the  collect 

There   is   one  dimensional   heat    flow   through    covers. 
There   is   a   negligible   temperature   drop   through   a   cover. 
There   is   one-dimensional  heat    flow    through  back  insulation, 
(body   for 

Temperature   gradients   around    tubes   can  be  neglected. 

The   temperature  gradients    in    the   direction  of    flow   and  between   the 

the   collector   absorbing  plate   is   negligible, 
plate   temperature  was   equal    to    the  mean  fluid    temperature 

izontal.      It  was   designed    to  be  able    to  drain  water   to  avoid 
ng.      It  was    "operated"  whenever  Tf  q,    the   collector   outlet 
,    was   less    than   210°F  and  was   greater    than  T      .,    the  collec- 


inlet   fluid    t 
the   collector 


The  major   design  paramete 
summed   up   in   Table   3.1. 


1    the   simulation 


i.l  Design  Parameters  of  the  Solar  Collectoi 


A       400; 600; 800; or  1000  ft   several  collector  areas  were  simulated 

(mc  )    10  Btu  hr"1"?"1  per  square  foot  collector  area 

(to)     0.81  assumed  independent  of  incident  angle 


neglect  bond  resistance,  see  Figure  3.2 

inside  tube  diameter 

forced  circulation 

thermal  conductivity  of  Aluminum 


1.64  x  10"J  ft  fin  thickness 

0.95  emittance  of  plat 


isulation  heat  loss  coeffici 


ENERGY  STORAGE 

Collected  solar  energy,  in  the  form  of  heated  water,  was  stored  in 
the  storage  tank.  The  stored  heated  water  was  supposed  to  supply  all  the 
heat  needed  for  house  and  hot  water  heating,  but  owing  to  the  shortage  of 
winter  sunshine  and  the  limitation  of  storage  tank  capacity,  a  supplementar 


The 


signed  to  operate  wit 


the  domestic  water  heating  circulation  loop  and  the  house  heating  circulation 
loop  (as  shown  in  Figure  3.3).   In  the  water  heating  circulation  loop,  the 
water  from  section  1  of  the  storage  tank  was  circulated  through  a  heat 
exchanger,  in  which  it  exchanged  heat  with  incoming  cold  water,  and  then 
returned  to  the  storage  tank  section  3.   The  temperature  of  the  domestic 
hot  water  leaving  the  heat  exchanger  was  assumed  10°F  lower  than  the  temp- 
water  was  not  as  high  as  the  minimum  acceptable  hot  water  temperature  (140°F) , 
the  hot  water  heater  supplied  energy  to  heat  it  to  that  temperature. 

The  hot  water  needed  was  assumed  to  be  uniform  at  47.3  lb/hr  from  7  AM 
to  7  PM  every  day.   The  temperature  of  the  domestic  cold  water  was  assumed 

was  lower  than  70°F,  the  circulation  to  the  heat  exchanger  was  stopped  and 

When  there  is  no  need  of  hot  water,  the  circulation  is  also  stopped.   When 
operated,  the  water  circulation  rate  in  this  loop  was  assumed  constant  value 
at  1000  lb/hr. 

In  the  house  heating  circulation  loop,  the  water  from  section  1  of  the 
storage  tank  was  circulated  through  coils,  heated  room  air  for  the  house, 
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was   designed   to   operate  at   a  minimum  temperature  of   120°F.      Once   the   temper- 
ature  of  water  returning    from   the   load  was  below   120°F,    the   controller  would 
stop   the   circulation  pump  and   the   supplementary    electric   furnace  would 
supply   all    the  house  heating   load.      The  water   circulation  rate,  when 
operated,    was    7000   lb/hr. 

A   three-section  water    storage   tank  model  was   used    in   the  simulations. 

fully  mixed.      The   energy  balance   equations    for    the   three   sections    can  be 

(m  c  )        !k 


(Ts,2        Ts,l)    +  F3(l°  Cp)LW(Ts,2        Ts,l) 


-    (UA)         (T  -   T   ) 


=  Fn(m   c    )    (T      ,    -   T 


s,2J  2ym  cp'nh^   s,3       s,2; 

"  V™  CpVTs,3  "  Ts,2> 

"    (UA)S;2(Ts,2  "   V  (12) 

c(Ts,2   "   Ts,3)   +  V™   Cp)HL(TLr"Ts,3) 
■    (UA>      ,<T      ,  -   V  (13) 


Table  3.2   Design  parameters  of  the  Storage  Tank 


4.17  Btu  F 

quare  foe 

4.17  Btu  F" 

same 

4.17  Btu  F~ 

same 

7000  Btu  hr"  °I 

stant 

1000  Btu  hr"  c 

stant 

11.6  Btu  hr"  ' 

a^con 

the  vali 
stant  foi 

8.49  Btu  hr"  * 

same 

8.49  Btu  hr"  " 

40''  F 

stant 

foot  of  collector  a 


"'    L"e    "°UVB    -4U.LXU,,,     ,Lr    *uu    ,LW    »u, 

tions.      Values   of   F   ,    F   ,    and   F3  are  determined  by: 

,1  when  T,   „    >   T,    .    and   I.  „   <   210°F 


when   T.       >   120  °F 


J  and   T  >   70°I 


c  Liier 


The    tank  was 


3  be  a  light  gauge,  vertical,  well-insulated, 
galvanized,  steel  cylinder  located  in  the  basement.  The  size  of  the  thei 
storage  tank  is  primarily  an  economic  decision.  Based  on  the  economic  st 
of   Tybout  and  Lof    [10],    the    thermal   equivalent  of   12.5   pounds   of   stored  w 


per   square   foot  of    collector 
are  summed   up   in  Table   3.2. 


Other  major   design  p 


DOMESTIC  HOT  WATER  LOAD 

The  hot  water  heating   load,    Q  was  assumed   constant   at   3784   Btu/h 


a  7  AM   to   7  PM   .      Since    the 


.0°F  lower    than    the    temperature   of  water   in  section   1 


of    the  storage   tank,    the  amount  of  \ 


:  heating   load   supplied  by   t 


"•3  j^)(l  ~)  (140-60)  (°F)    =   3784  ^        (7  AM   to   7   PM) 

(16A) 


QBS  -  (47.3)(l)(Thot  -  60) 
where: 


And,  the  temperature  of  water  returning  from  the 
to  section  3  of  the  storage  tank  was  calculated  b 


HOUSE  HEATING  LOAD 

from  the  building.   The  heat  losses  may  be  divided  into  two  groups:   (1)  the 
transmission  losses,  or  heat  transmitted  through  the  confining  walls,  floor, 
ceiling,  glass  or  other  surfaces;  and  (2)  the  infiltration  losses,  or  heat 
required  to  warm  outdoor  air  which  leaks  in  through  cracks  and  crevices, 
around  doors  and  windows,  or  heat  required  to  warm  outdoor  air  used  for 
ventilation. 

In  this  study,  a  simple  space  heating  model,  the  energy  per  degree-time 
model  [12],  was  used  in  estimating  the  heating  load.  The  house  heating  load 
QHL  was  expressed  as  a  linear  function  of  the  difference  between  the  inside 
temperature  of  the  building  and  the  ambient  temperature, 

QHL  -  (U)(A)(T.  -  Ta)  (20) 

The  design  indoor  temperature  was  70°F.   The  value  of  uA  for  the  Simula- 
Table  3.3  and  was  equal  to  1014  Btu  hr"1"?"1. 
(1)   The  residence  was  assumed  to  be  heatf 
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(2)  The  basement  was  assumed  to  need  no  heating  or  cooling. 

(3)  The  transient  heat  from  the  storage  tank  warms  the  air  near  the  base- 
ment ceiling  sufficiently  to  make  it  unnecessary  to  make  an  allowance 
for  floor  heat  loss  from  rooms  located  over  the  basement. 

(A)   The  humidity  difference  between  the  indoor  air  and  the  outdoor  air 

(5)  The  heat  supplied  by  persons,  lights,  motors,  machinery,  and  sunshine 
were  negligible. 

(6)  The  wind  speed  was  assumed  to  be  a  constant,  15  M.P.H. 

Once  the  heating  load  was  calculated  by  equation  (20),  the  determina- 

v  p'hl 

HOUSE  COOLING  LOAD 

by  heat  flow  through  structural  components  and  by  air  leakage  or  ventila- 
tion.  In  this  study,  the  calculation  of  cooling  load  was  based  on  the 
estimation  methods  of  ASHRAE  [30]  and  Allen  [31].   The  residence  was 
assumed  to  be  cooled  as  a  single  zone,  and  usually  conditioned,  for  24 
hours  per  day,  every  day  of  the  cooling  season. 

The  total  residential  cooling  load  is  the  sum  of  the  sensible  heat  load 
and  the  latent  heat  load.  A  sensible  heat  load  is  considered  to  occur  when 
there  is  a  direct  addition  of  heat  to  the  enclosure  by  any  one  or  all  of 

load  is  considered  to  occur  when  there  is  an  addition  of  water  vapor  to  the 


the  structure,  the  occupancy,  and  the  equipment;  the  latent  portion  of  a 

of  the  calculated  sensible  heat  load.   The  calculation  of  residential  cooling 
load  is  therefore  invariably  the  calculation  of  the  sensible  heat  load. 
In  calculation  of  the  sensible  heat  load,  the  following  items  must 
be  considered: 

(1)  heat  gain  through  walls  and  roof 

(2)  heat  gain  through  windows 

(3)  heat  gain  due  to  infiltration  and  ventilation 

(4)  heat  gain  due  to  occupancy  loads 

The  detailed  sensible  heat  load  calculations  of  the  above  four  items 
are  as  follows: 

(1)   heat  gain  through  walls  and  roof 

The  heat  gain  through  a  wall  may  be  expressed  by  the  equations  proposed 
by  Faust,  Levine,  and  Urban  [32], 

"wall  =  Ht,„all  +  Hr,wall  (22) 

Ht,wall  "  (A„all)(U„all)(Ta  "  V  (23) 

"r.wall  =  (A„ii»(P)(')(I.)  (2M 

The  values  of  P,  the  portion  of  the  absorbed  solar  radiation  that  is 

P  -  (0.23)(l 


.-a  1 1 


and  Eq.  (24)  becomes 

Hr,wall  =  (0-23)(Awall)(Uwall)(a)(Ia)  (26) 

The  absorption  coefficient,  a,  was  assumed  to  be  0.4  for  the  light  colored 
wall  surfaces.  Values  of  Ifl,  the  solar  radiation  intensity,  vary  with  time 
and  the  direction  the  walls  face. 

The  calculation  of  heat  gain  through  the  roof  is  almost  the  same  as 
the  calculation  of  heat  gain  through  a  wall  except  that  it  was  assumed 
that  the  collector  was  located  on  the  roof,  and  thus  shaded  the  roof, 
eliminating  heat  gain  by  solar  radiation.   The  equation  for  calculating  the 
heat  gain  through  the  roof  was: 

Hroof  =  (Aroof)(Uroof)(Ta  "  V  (27) 

According  to  Allen  [31],  the  thermal  capacity  of  the  outer  walls  and 
roof  has  a  marked  effect  by  causing  a  time  lag  in  the  penetration  of  heat 
to  the  room.  This  lag  may  be  as  much  as  3  hour  for  a  6-in.  concrete  wall 
and  10  hour  for  a  22-in.  brick  and  tile  wall.   So  that,  in  applying  the 

It  was  assumed  to  be  4  hours.   The  heat  gain  through  walls  and  roof  were 
calculated  as  in  Table  3.4. 

The  heat  gain  through  windows  may  be  expressed  by  a  simple  equation: 

"w  =  AwFwIa  +  AwUw(Ta  -  Tj)  (28) 

The  solar  heat  gain  coefficient  Fw  is  a  characteristic  of  each  type  of 

tance  of  the  glazing  material  are  dependent  upon  the  incident  angle.   But 
in  this  study,  it  was  assumed  that  transmittance  and  absorptance  of  the 
glazing  material  were  constant,  so  that,  the  value  of  Fw  was  also  assumed 
independent  of  the  incident  angle. 


The  calculation  of  Fw  was  based  on  the  Shading  Coefficient  described 
n  ASHRAE  [30].   The  Shading  Coefficient,  S.C.,  is  represented  by  the 


S.C.  =  ■ 


Dlar  Heat  Gain  of  Fenestration 


Solar  Heat  Gain  of  Double-Strength  Glass  W 

In  terms  of  the  Solar  Heat  Gain  Coefficient,  F  ,  the  Shading  Coeffic 
:  ratio  of  Fw  for  the  fenestration  and  F  for  D.S.A.  glass  which,  for 


F  of  D.S.A.  =0.87 


(30) 


For  insulating  glass  with  indoor  shading  by  opaque,  white  roller  shades, 
the  Shading  Coefficient  was  determined  to  be  0.25.   Based  on  this  value,  the 
Solar  Heat  Gain  Coefficient  Fw  of  the  window  was  calculated  by  Eq.  (30)  to 
be  0.22.   The  Heat  Gain  Through  Windows  was  calculated  as  in  Table  3.5. 

Natural  air  leakage  in  residential  structures  is  much  smaller  in  summer 
than  in  winter.  According  to  ASHRAE  [30],  Infiltration  and  Ventilation  air 
flow  is  usually  assumed  at  a  rate  of  about  one  air  change  per  hour  in  summer. 

The  calculation  of  heat  gain  due  to  infiltration  and  ventilation  was 
based  on  the  equation  (reference  (29)): 

Hj.  =  (0.018)  (VQ)(Ca)(Ta  -  T±)  (31) 

Since  VQ  =  18000  ft3,  Ca  =  1,  and  T.  =  75°F,  the  heat  gain  due  to  infiltra- 
tion and  ventilation  was  calculated  as  Hj  =  (324)(Ta  -  75). 

Even  though  occupant  density  is  usually  lower  in  residences  than  in 


usually  assumed  to  be  225  Btu  per  hour  of  sensible  heat.   Appliance  load 
can  be  limited  to  the  kitchen  in  most  cases  and  is  assumed  to  be  1200  Btu 
per  hour  of  sensible  heat.   In  the  simulation  model,  the  number  of  occupants 
was  estimated  to  be  6  and  the  appliance  loads  were  assumed  to  be  at  5-6  AM, 
11-12  AM,  and  5-6  PM. 

load  calculation,  total  load  is  usually  calculated  as  1 . 3  times  the  calculate 
sensible  load. 

WEATHER  DATA 

The  weather  data  needed  in  the  simulation  study  were  hourly  solar  radia- 

calculating  collector  heat  loss  coefficient).   These  weather  data  were 
provided  by  the  Department  of  Physics,  Kansas  State  University. 

The  hourly  solar  radiation  data  provided  were  the  direct  and  diffuse 
solar  radiation  in  gm-cal  per  cm2  of  horizontal  surface.   The  actual  hourly 
solar  radiation  impinging  on  the  tilted  solar  collector,  the  vertical  walls 
and  windows  were  calculated  by  first  determining  the  solar  angle  for  every 

According  to  ASHRAE  [30],  the  incident  angle  6  is  related  to  the  solar 
attitude  B,  the  wall-solar  azimuth  y,    and  the  tilt  angle  e,  by: 


ns  needed  for  the  calculation  of  cos6  are: 

B  =  cosL  cos6  cosH  +  slnL  sin6  (3: 

t   =   cos6  sinH  /  cosB  (3< 

0.25  x  (Number  of  minutes  from  local  solar  noon)         (31 
r  declination  S,  which  is  a  function  of  the  date  (see  Table 
jmed  to  be  a  constant  for  each  month.   The  relationship  among 
Les  is  as  shown  in  Figure  3.4. 

5  to  Eq.  (30),  when  a  surface  is  horizontal,  e=0°,  then 
3H  -  sinS  (36 

rtical  surface,  e=90° 


Also,    foi 
(e-49"), 


the   collector   facing  south    ( 


Jith   a    tilted   angle  of  i 


i   solar  radiation   impinging  on   the   collector   can  be   calculated  by 
wing   equation: 


The  value 
facing   vet 


)f   y,    the  wall-solar  azimuth,    in    the  above   equation    for  differently 
:ical  walls   and  windows  are  given   in  Table   3.7. 


Table  3.6  Values 


f  Solar  Declinatic 


DECLINATION  ,  degre 


Aug  2: 
Sep  2] 


SOLAR  ALTITUDE  ,  ^3  =  L  QOH 
ZENITH  ANGLE  ,  Z  -  Z.  QOV 
HALL  AZIMUTH ,       Us    =  /_  SOP 


SOLAR  AZIMUTH,  <f>  =  Z  SOH 
INCIDENT  ANGLE,  9V=  £QOP 
WALL-SOLAR  AZIMUTH,     \  =  Z.  HOP 


Figure  3.4     Definition  of   Solar  Angles 


l   for  Different  Facing  V 


DIRECTION 

r,de5^ 

NOTES 

(h     -  90  AM 
(p    +   90  PM 

2.  Treat 
r  as 

3.  If  f i 

alculated  by  Eq.(34). 

s  greater  than  90* 
rface  is  in  the  shade. 

».. 

<£  +  90  AM 
<p    -   90  PM 

— 

<£ 

— 

greater  than  90' 

The  derivations  of  Eq.  (39)  and  Eq.  (40)  were  based  on  the  assumption 
that  most  of  the  diffuse  radiation  comes  from  an  apparent  origin  near  the 

Under  this  assumption,  the  horizontal  radiation  was  treated  as  though  it 
was  all  beam  radiation.  Although  this  approximation  is  best  suitable  on 
clear  days,  it  should  not  cause  large  errors  in  the  simulations. 

AUXILIARY  ENERGY 

In  this  model,  the  auxiliary  energy  for  house  heating  and  for  hot 
water  heating  was  supplied  separately  by  an  electric  furnace  and  an  electric 
hot  water  heater.   The  furnace  was  assumed  to  be  able  to  supply  all  the  heat 
load  without  having  any  time  lag  once  the  controller  stops  the  circulation 
of  heated  water  from  the  storage  tank  to  the  load  (when  TLr  <  120°F) .   A 
similar  assumption  was  made  on  the  hot  water  heater  -  i.e.,  that  it  could 
supply  energy  to  supplement  the  solar  energy  to  meet  the  hot  water  load  at 

and  for  hot  water  heating  was  estimated  in  the  simulation  by  the  following 
equations: 

(Qm  when  I   <  120  °F 


-  3784  when  T.^  <  70°F  and  hot  water  is 

(47.3)(140-T   )   when  60°F  <  T    <  140°F  and  hot  s 


calculated   by: 

<W  =  Z((W(At)  w 

^AUW^W^  C44) 

where  At    is    the    time   interval;    one  hour  was    chosen   in   the   simulations. 

In   the   calculations   above,    the   electric    energy   consumed  by    the   circula- 
tion pumps  was   neglected.      The   same  assumption  was  also  made   in  MODEL  B 
and  MODEL  C. 

of    the  air  conditioner  was    estimated   by    the  value  of   CO. P.    (Coefficient 
of   Performance).      This  will  be  discussed    in    the   section  -   HEAT  PUMP.      (It 
was   assumed   in    this    study    that    the   CO. P.    of   an  air   conditioner  and    a  heat 
pump  was    the  same  when    they  are  operated    in   the   same  weather    conditions.) 

MODEL  B   -    SOLAR-ASSISTED  HEAT  PUMP   FOR   HEATING  AND   COOLING 

A  schematic   diagram  of    the  heating   portion  of    this  model    is   shown   in 
Figure   3.5.      Unlike   in  MODEL  A,    a  heat   pump  was    installed   for  both   heating 

the  necessary  pumps,    fan,    and   controllers  were  used   for  house  heating.      The 

It  was   assumed    to  be  just    like   an  air   conditioner  when  used    for   cooling. 

Except   for    the  heat   pump,    the   energy   storage,    and    the  auxiliary   energy, 
all   component    simulations  were   the   same   for    this   model   as  was    the   case 
for    the   previous  model.      The   energy   storage,    the  heat   pump  and    the  auxiliary 
energy  were  modeled   as   follows. 


Figure  3.5  Solar-A 


:  Pump  Heating  Syst 


ENERGY  STORAGE 

In  this  model,  the  water  in  the  storage  tank  was  again  designed  to 
operate  in  two  major  circulation  loops:   the  hot  water  heating  loop  and 

as  in  the  previous  model,  and  will  not  be  discussed  again.   In  the  house 
heating  circulation  loop,  there  are  two  routes  for  water  circulation. 

heating  and  then  returned  to  the  storage  tank  section  3.   The  second  route 
is:   water  from  the  storage  tank  section  1  was  circulated  through  the  heat 
pump  suction  side  heat  exchanger  in  which  it  heated  the  refrigerant  and 
then  returned  to  the  storage  tank  section  3. 

When  the  temperature  of  water  returning  from  the  load  (heating  coils) , 
TLr,  is  higher  than  120°F,  the  circulation  pump  will  circulate  the  water 
through  the  first  route,  and  the  stored  heated  water  is  responsible  for 
all  the  house  heating  load.   When  TLr  is  lower  than  120°F  and  the  temperature 
of  water  returning  from  the  heat  pump  heat  exchanger,  TLh,  is  higher  than 
40°F,  the  circulation  pump  will  circulate  the  water  through  the  second 
route.   In  this  case,  the  stored  heated  water  is  used  as  a  heat  source  for 
the  heat  pump,  and  the  heat  pump  is  responsible  for  all  or  part  of  the  house 
heating  load.   When  TLr  <  120°F  and  TLh  <  40°F,  the  controller  will  stop 
the  circulation  pump,  and  the  electric  furnace  will  supply  all  the  heat  for 

Such  operation  can  be  expressed  by  previously  derived  equations,  Eq.  (11) 
(12) (13) (14)  and  (16),  except  that  Eq .  (15)  is  replaced  by: 


when  TLr  >  120°F 

when  TLr  <  120°F  and  TLh  >  40°F  (with  T^      („ 
replaced  by  TLh  in  Eq.  (13)) 
when  TLr  <  120°F  and  TLh  <  40 °F 
Hie  calculation  of  T^  was  by  Eq .  (21),  and  the  calculation  of  TLh  will  be 
discussed  in  the  next  section.   Other  portions  were  modeled  the  same  as  in 
MODEL  A. 

HEAT  PUMP 

TLh  >   40°F>  and  was  °Perated  in  the  cooling  season  just  like  a  conventional 
air  conditioner.   A  heat  pump  heating  cycle  is  shown  in  Figure  3.6.   The 

The  heat  pump  capacity  was  assumed  to  be  a  function  of  the  heat  pump 
CO. P.  (Coefficient  of  Performance).   The  ideal  heat  pump  CO. P.  is  defined 

C-°-P-I  =  work  required  =  T^T  <«> 

According  to  Kemler  and  Oglesby  [33],  the  ratio  of  the  CO. P.  of  a 
commercial  compressor  to  the  ideal  CO. P.  is  just  about  one  half.   This 
ratio  is  a  weak  function  of  both  the  suction  and  the  condensing  temperature 
In  this  heat  pump  model,  a  constant  of  0.53  was  assumed  to  be  the  ratio  of 
a  commercial  heat  pump  CO. P.  compared  to  the  ideal  CO. P.   Then,  the  heat 
pump  CO. P.  was  calculated  by: 


n^" 


£4= 


Figure  3.6   The  Heat  Pump  Heating  Cycle. 
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Figure  3.7    The  Heat  Pump  Cooling  Cycle 


heat  rejected  =  heat  abosrbed  +  net  work  input  (48) 

finally  becomes  the  output  heat  of  the  heat  pump.   By  assuming  that  all 

pressor  into  work,  the  auxiliary  energy  required  (electric  energy  consumed 
by  the  compressor)  is  thus  equivalent  to  the  heat  added  by  the  compressor 
(or  work  input) . 

capacity.   It  was  determined  by  rearranging  Eq.  (47), 

(the  heat  pump  capacity,  cy  =  (C .0. P . ) (work  required) 

=  (CO. P.)  (horsepower  input)     (49) 
The  horsepower  input  was  treated  as  a  constant  which  depended  on  the  size  of 
the  heat  pump  compressor.   Once  the  compressor  size  was  fixed,  the  heat 
pump  capacity  was  therefore  a  function  of  CO. P.  only.   In  this  heat  pump 
model,  the  suction  temperature  ^  was  assumed  to  be  10°F  below  that  of  the 
temperature  of  water  in  the  storage  tank  section  1.   The  condensing  tempera- 
ture was  designed  to  be  120°F. 

When  the  heat  pump  was  operated  for  house  heating,  the  procedures  of 

QAUX'  iTlstantaneous  hourly  auxiliary  energy  required,  are  as  follows: 

(1)  Calculate  CO. P.  by  Eq.  (47) 

(2)  Determine  Ch  by  Eq.  (49) 

(3)  If  Ch  >  QHL,  this  means  that  the  operation  of  the  heat  pump  alone 


(4)   If  Ch  <  QHL>  the  electric  furnace  is  operated  to  supplement  the 


XAUX   CO. P.    VVHL 


operate  in  the  same  manner  as   a   conventional  air  conditioner.   A  heat  pump 
cooling  cycle  is  shown  in  Figure  3.7.   In  this  case,  the  condensing  tempera- 
ture was  designed  to  be  55°F.   The  value  of  the  CO. P.  was  calculated  by 
Eq.  (47),  and  the  hourly  auxiliary  energy  required  (electric  energy  consumed 
by  the  compressor)  was  calculated  by: 

house  cooling  load    when  cooling  is  required 
Q    .  I      C.O.P.  -  1         (cooling  load  >  0) 

I  0  otherwise  (54) 

(This  assumed  that  the  size  of  the  heat  pump  or  the  air  conditioner  was 
large  enough  to  meet  all  the  house  cooling  load.) 

AUXILIARY  ENERGY 

water  heating  in  this  simulation  model: 

(1)   Electric  energy  required  to  operate  the  heat  pump  compressor. 


required  for  house  heating  were  as  follows: 

(1)  When  TLr  >  120°F,  all  the  house  heating  load  was  supplied  by  the 
stored  heated  water,  and 

<W  "  °  (55) 

(2)  When  TLr  <  120°F,  and  TLh  >  40°F,  QAUX  is  expressed  as  in  Eq .  (50) 
and  (52). 

(3)  When  TLr  <  120  °F  and  T   <  40  °F, 

w  ■  V  (56) 

(4)  The  cumulative  auxiliary  energy  required  for  bouse  heating  was 
calculated  by  adding  all  the  hourly  auxiliary  energy  requirement 
terms  during  the  simulation, 

WZ((W(At)  (57) 

The  hourly   auxiliary   energy   required    for   house   cooling  and    for  hot 
water  heating  was   calculated  by   Eq .    (54)    and  by   Eq .    (42).      The   long-term 
cumulative  auxiliary   energy  required   for    these   two   items  was  again  calcu- 
lated by  adding   all    the  hourly    terms   during   the  simulation. 


chill    the  water    in    the   storage   tank  during    the  summer  night,    and    this 


The  heat  pump  was  designed  to  start  operating  at  midnight  every  night 
until  the  temperature  of  water  in  the  storage  tank  section  1  reached  40°F. 
The  chilled  water  in  the  storage  tank  section  3  was  circulated  through 
cooling  coils  to  meet  the  cooling  load  once  the  house  cooling  was  needed, 


and  then 


)  the  storage 


to  chill  the  water  so  that  the 
The  energy  balance  equations  fo: 


For 


i  1) 


"Vs.l        dt  r5        p- 


Lon  3! 


V'"V 


of 


ink  section  1.      When   the   tei 

lerature  of  water    for   cooling  would  not 
.ng   temperature  of   55°F. 
ooling  operation   is   shown   in  Figure   3.! 
;e  storage   tank  were  as    follows: 

CL(TCL  *  ^.l5   +VS)CH(Ts,2-Ts>l) 
..l'Wl'  (5i 

(T  -  T         )    +   F    fmr    1        fT  -  T         ■> 

CLUs,l        s,2;        VmVcHUs,3      rs,2J 

Sj2(TB-Ts>2)  (5S 

p'CL^s^"^5    +  F6(mcp)CH(TCH"Ts,3) 


i 


>oling  was   needed 


F      = 


I  H     "  I 


'a-S.3'  (S)CL 


In  the  equation  above,  the  value  of  C  and  the  value  of  CO. P.  were  determ 
by  Eq.  (49)  and  by  Eq .  (47) . 

Finally,  the  hourly  auxiliary  energy  required  (electric  energy  consu 
by  the  heat  pump  compressor)  was  calculated  by: 


C  .  0  .  P  . 


len  the 


temperature  was  assumed  to  be  10°F  below  the  temperature  of  water  in  the 
storage  tank  section  1.   The  basement  temperature  was  assumed  to  be  70°F. 
The  values  of  (mc  )    and  (mc  )CH  were  designed  as  7000  and  12000  Btu/hr°F. 
The  simulation  results  of  the  three  models  are  reported  in  CHAPTER  5. 


CHAPTER  4.   CONTINUOUS  SYSTEM  MODELING  PROGRAM 

role  in  system  simulation  and  in  the  solution  of  nontrivial  mathematical 
equations.   Computer  languages  such  as  FORTRAN,  BASIC,  and  PL/1  have  pro- 
vided users  with  a  convenient  method  for  communic sting  with  the  computer  on 
a  wide  range  of  problems.   In  man}'  cases,  specialized  computer  programs 

problems.   CONTINUOUS  SYSTEM  MODELING  PROGRAM  (CSMP)  is  one  of  those  kind 
of  programs.   It  is  most  suitable  for  the  simulation  of  physical  systems. 

The  CONTINUOUS  SYSTEM  MODELING  PROGRAM  (developed  by  the  IBM  Company) 
was  designed  to  allow  users  to  simulate  all  types  of  physical  systems  with 
a  minimum  of  programming  difficulty.  It  is  an  application-oriented  program 
in  that  it  is  specially  written  for  scientists,  engineers,  and  analysts  who 
are  involved  in  work  that  requires  the  solution  of  ordinary  differential 
equations  or  in  simulating  a  system  that  has  been  modeled  as  a  block 
diagram. 

In  developing  the  program,  emphasis  was  placed  on  simplified  input  data 

directly  describe  the  mathematical  equations  or  physical  variables  of  the 
problem.   The  program  which  is  an  extension  of  FORTRAN  language  includes 
a  basic  set  of  functional  blocks  which  can  represent  the  components  of  a 
continuous  system  and  accepts  application-oriented  statements  defining  the 
connections  between  these  functional  blocks.   Included  in  the  basic  set  are 

many  special  purpose  functions  like  delay  time,  zero-order  hold,  dead  space, 
limiter  functions,  and  so  on  Those  functions  really  make  the  CSMP  an 


54 
attractive  language  in  simulating  complex  physical  systems.   For  example, 

involved  in  specifying  either  the  type  or  any  of  the  details  of  the  numerical 
integrating  method.   This  simplifies  the  programming  and  reduces  the  chances 
of  making  mistakes.   In  effect,  CSMP  allows  the  user  to  concentrate  on  the 
details  of  the  physical  system  rather  than  the  usual  concerns  of  numerical 
analysis  and  programming. 

Using  the  CSMP  program  is  similar  in  many  respects  to  using  the  elec- 
tronic analog  computer.   Many  of  the  special  CSMP  statements  perform  the 
same  function  as  typical  analog  computer  components.   In  fact,  CSMP  digital 
simulation  has  many  advantages  over  the  use  of  an  analog  machine  in  that 

more  accurate;  (3)  is  easier  to  program;  and  (4)  has  much  greater  capability 

In  addition  to  simulating  solar  processes,  applications  in  which  CSMP 

parameter  estimation,  studies  of  blood  circulation  and  other  physiological 
processes,  studies  of  chemical  refineries,  natural  gas  transmission,  process 
control,  investigation  of  aircraft  landing  and  take-off,  plant  growth, 
natural  resources  management  and  other  industrial  dynamics.   Anyone  working 
in  a  field  of  science  who  has  a  need  to  simulate  a  system  or  solve  ordinary 
differential  equations  will  quickly  recognize  the  power  of  CSMP. 

CSMP  is  most  useful  for  small  and  medium-size  simulations.   For  extremely 
large  and  complex  problems,  the  simulation  may  be  limited  by  program-size 
restrictions.   Also,  because  of  the  flexibility  in  programming,  CSMP  simula- 
tions may  require  slightly  more  computer  time  than  custom-developed  programs. 


CHAPTER  5.   S' 


HEATING  RESULTS 

energy  quantities.   These  include  the  cumulative  hot  water  heating  load, 
Q^;  the  cumulative  house  heating  load,  QHL;  the  cumulative  auxiliary  energy 

Two  simulation  programs  were  written  using  CSMP  (Continuous  System 
Modeling  Program)  to  calculate  the  above  mentioned  terms.   The  programs 
were  written  according  to  the  component  models  described  in  Chapter  3. 
The  equations  among  the  component  models  were  coupled  and  were  well  arranged 
so  that  they  gave  the  correct  simultaneous  solution.   The  simulation  algor- 


jis  were  as  follows: 

(2)  Read  the  weather  da 

(3)  Calculate  the  solar 
Eq.  (39). 

(4)  Calculate 

(5)  Determine 


is,  such  as  initial  tank  temperature, 


(7)   Calculi 

*WL  and  %L 


illector  heat  loss  coefficient,  U  ,  by  Eq.  (6). 
:he  collector  outlet  fluid  temperature,  Tf   ,  by  Eq .  (1) 
lat  the  collector  inlet  fluid  temperature  is  the  same 
on  3  fluid  temperature  of  the  storage  tank.) 

ng  load,  QHL,  by  Eq .  (16A)  and  by  Eq .  (20). 


(10)  Ci 

(11)  Ca: 


(14)  D( 

(15)  Ca: 


xhanger,  1^,  by  Eq.  (19). 

ite  the  temperature  of  water  returning  from  the  house 

;  load,  TLr,  by  Eq.  (21). 

ite  CO. P.    of    the   heat   pump  by  Eq .    (47). 

'   Eq.    (51). 

te  hourly   auxiliary   energy   required    (QAUW  and    Q^)    by 
1)    and    (42)    for  MODEL  A   or  by  Eqs.    (42),    (50),    and    (52) 
MODEL  B  and   C. 

te  cumulative  auxiliary  energy  required,  \m  and  QAUX, 
(43)    and  by  Eq .    (44). 

ne  F1,  F2,  and  F3  by  Eq .  (14),  Eq .  (15),  and  Eq.  (16). 
te  the  temperature  of  water  in  sections  1,  2,  and  3  of 
rage   tank   (T^,    T,^,    and   T     3)    by  Eqs.    (11),    (12),    and 


(13) 

(16)  Calculate  collector  plate  temperature,  T  ,  by  Eq .  (7). 

(17)  Repeat  steps  2-17  for  the  next  time  interval. 

Based  on  the  equations  derived  in  Chapter  3,  the  calculation  results  of 
the  cumulative  auxiliary  energy  requirements  have  units  of  Btu.   Since  all 

was  assumed  that  the  conversion  efficiency  from  electricity  into  heat  was 
1.0,  the  termal  energy  provided  by  the  auxiliary  systems 


f   energy 


equired)  was  thus  equivalent  to  the  elect 


ergy  c 


MODEL  B  and  C  only. 


Since  the  simulation  of  MODEL  B  required  the  value  of  the  heat  pump 

heat  pump  needed  to  be  chosen  before  the  simulation  was  made.  Based  on 
calculations  of  the  cooling  load  for  the  hottest  month  (see  Figure  5.1), 

large  enough  to  meet  the  cooling  requirements  of  the  house  in  the  summer, 
therefore,  a  3-ton  air  conditioner  was  used  in  MODEL  A*  and  a  3-ton  heat 
pump   was  used  in  MODEL  B  '  and  MODEL  C  . 

In  order  to  minimize  the  computer  costs,  only  the  first  seven  days  of 
each  month  requiring  heating  (Oct.  -  Apr.)  were  simulated.   The  weather 
data  of  the  first  seven  days  of  each  month  were  punched  in  the  computer 

simulation  was  made  by  using  the  weather  data  at  the  end  of  one  seven  day 
period  as  the  initial  values  for  the  beginning  of  the  next  seven  day  period. 
The  performance  of  this  seven  day  simulation  was  then  multiplied  by  a  factor 
(31/7  for  October,  December,  January,  and  March;  30/7  for  November  and  April; 
and  28/7  for  February)  to  represent  a  whole  month's  result. 

The  computed  monthly  and  yearly  system  performance  for  heating  for 
MODEL  A  and  MODEL  B  with  collectors  of  400,  600,  800,  and  1000  ft2  are  shown 
in  Tables  5.1  through  5.4  (MODEL  C  has  the  same  performance  as  MODEL  B.) 

supplied  by  these  collectors  for  both  models 


3ad  (0-AW  +   W  f"C 


i   against 


sa  in  Figure  5.2.   The  figure  shows  clearly  that  MODEL  B 
L  performance  than  MODEL  A. 
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Some  dynamic  results  during  the  simulation  are  also  plotted  in  Figures 
5.3  through  5.28.   In  these  figures,  the  units  of  time  are  days.   Since  only 
the  first  seven  days  of  each  month  were  simulated,  the  performance  between 
0  and  7  days  stood  for  the  performance  for  October,  and  the  performance 
between  7  and  14  days  in  these  figures  stood  for  the  performance  for 
November,  etc. 

Figures  5.3  through  5.10  show  the  hourly  auxiliary  energy  required  for 
hot  water  heating  for  each  model  with  collector  areas  of  400,  600,  800,  and 
1000  ft2  during  the  simulation.   Every  figure  shows  that  more  auxiliary 

and  the  least  amount  of  auxiliary  energy  was  needed  in  November  (between 
days  7  and  14).   By  comparing  these  figures,  we  can  see  that  as  the  collector 


ergy  requii 


heating  for  MODEL  A  is  always  less  than  the  amount  required  for  MODEL  B 

in  the  storage  tank  for  MODEL  A  is  most  often  higher  than  that  for  MODEL  B 
during  the  simulation  (as  shown  in  Figures  5.27  and  5.28).  The  figures  for 
MODEL  B  (Figures  5.7  through  5.10)  also  show  that  domestic  hot  water  load 
is  often  supplied  solely  by  the  auxiliary  heater  because  the  auxiliary  ener 
value  shown  in  the  figures  often  reaches  its  maximum  (3784  Btu/hr) .  This 
indicates  that  the  storage  tank  temperature  is  often  below  70°F  for  MODEL  B 
during  the  simulation. 

Figures  5.11  through  5.18  show  the  hourly  energy  required  for  space 
heating  for  each  model  with  each  of  the  four  collector  areas.   All  the 
figures  show  that  the  highest  hourly  auxiliary  energy  required  for  space 


heating  happens  at  day  20  (in  December).  However,  the  month  which  requires 
the  greatest  amount  of  auxiliary  energy  is  January  since  the  area  under  the 
curve  is  the  largest  in  this  month.   These  figures  also  show  that  the  amount 

By  comparing  the  figures  for  the  two  models,  MODEL  B  is  shown  to  need  less 
auxiliary  energy  than  MODEL  A  for  the  same  collector  area.  This  indicates 
that  the  solar-assisted  heat  pump  system  has  a  better  thermal  performance 
than  the  conventional  solar  energy  system. 

The  cumulative  total  auxiliary  energy  required  was  plotted  against 
time  in  Figures  5.19  through  5.26.   These  figures  show  clearly  that  the 

energy  required.  These  figure;  also  show  that  the  amount  of  the  cumula- 
tive auxiliary  energy  required  for  MODEL  B  is  always  less  than  that  for 
MODEL  A  for  the  same  collector  size  during  the  simulation. 

was  also  plotted  against  time  in  Figures  5.27  and  5.28.   These  figures  show 
that  the  water  temperature  in  the  storage  tank  varies  during  the  simulation 
between  100°F  and  210°F  for  MODEL  A  and  between  40°F  and  210°F  for  MODEL  B. 
The  weather  data  during  the  simulation  and  the  simulation  programs  for 
MODEL  A  and  MODEL  B  are  shown  in  the  Appendix. 

An  entire  month's  simulation  for  MODEL  A  in  January  was  also  made. 
The  attempt  is  to  compare  the  performance  of  this  month  simulation  with 

seven  day  simulation.   The  results  are  shown  in  Table  5.5.   It  shows  that 


Table  5.5  Comparison  Between  the  Simulated  Monthly  Performance  and  the 
Estimated  Monthly  Performance. 
(MODEL  A,  Collector  Area  ■  600  ft2,  January) 


Energy  Quantities  (  x  10'   Btu  ) 


0.0%         -0.8%         12.87.        17.5 

ntire  month's    simulation   results. 

Lrst   seven  day's    simulation   results   multiplied   by  31/7. 
on  simulated    results    =  100%. 
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COOLING  RESULTS 

The  cooling  result  needed  for  the  cost  analysis  is  the  cumulative  auxi- 
liary energy  required  (total  electric  energy  consumed  by  the  compressor) 

required  was  calculated  by  adding  all  the  instantaneous  hourly  auxiliary 
energy  required  values  which  were  estimated  by  the  values  of  the  CO. P. 
and  were  calculated  by  Eq.  (49)  for  MODEL  A  and  B  and  by  Eq.  (59)  for 
MODEL  C 

The  simulation  result  showed  that  the  cumulative  auxiliary  energy 
required  for  MODEL  A  and  B  was  less  than  that  for  MODEL  C  (see  Figure  5.29). 
This  result  was  caused  by  the  fact  that  the  CO. P.  value  for  MODEL  C  during 
the  simulation  was  mostly  lower  than  that  for  MODEL  A  and  MODEL  B  (see 
Figure  5.30).   This  means  that  the  operation  of  a  heat  pump  during  the 

MODEL  C  has  a  lower  suction  temperature  since  it  was  assumed  that  the 
suction  temperature  for  MODEL  C  was  10°F  below  the  section  1  temperature 
of  the  storage  tank.   The  storage  tank  temperature  for  MODEL  C  is  varied 
during  the  simulation  between  55°F  and  40°F,  so  that,  the  averaged  suction 
temperature  for  MODEL  C  is  about  37°F  instead  of  55°F  for  MODEL  A  and  B. 
This  means  that  MODEL  C  will  have  better  performance  only  when  the  temper- 
ature difference  between  daytime  and  nighttime  is  greater  than  20°F.   This 
is  not  the  case  for  Manhattan's  weather,  so  that  MODEL  C  showed  less 
efficiency.   Table  5.6  clearly  shows  this  fact. 
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In  Table  5.6,  the  averaged  condensing  temperatures  were  chosen  as 
explained  because  the  operation  of  the  heat  pump  in  MODEL  B  was  mainly 
during  the  daytime  from  9  AM  to  about  9  PM,  and  the  operation  of  the  heat 
pump  in  MODEL  C  was  mainly  from  midnight  to  about  7  AM.   It  is  evident  fron 
Table  5.5  that  MODEL  C  would  consume  more  auxiliary  energy  than  MODEL  B, 
therefore,  MODEL  C  will  not  be  discussed  further. 

In  the  next  chapter,  a  cost  analysis  is  made  to  compare  the  economic 
feasibility  of  MODEL  A  and  MODEL  B  with  a  conventional  heating  and  cooling 
system  which  uses  an  electric  furnace  for  house  heating,  an  electric  hot 
water  heater  for  hot  water  heating,  and  an  air  conditioner  for  house 

the  basis  (operation  and  initial  costs  of  a  conventional  system),  and  since 
energy  required  for  house  cooling  is  not  needed  and  was  not  calculated  in 


CHAPTER  6.   COST  ANALYSIS 

feasibility  of  the  two  solar  energy  systems.   The  method  of  analysis  was 
to  determine  the  annual  dollar  savings  of  the  systems.   The  term  "annual 

using  a  solar  system  for  heating  and  cooling  instead  of  using  a  conventional 
heating  and  cooling  system.   It  was  calculated  by  subtracting  the  annual 
costs  of  each  system  from  that  of  a  conventional  system.   The  conventional 
system  chosen  in  this  study  is  one  which  uses  an  electric  furnace  for  space 
heating,  an  electric  heater  for  domestic  hot  water,  and  a  3-ton  air  conditio 

The  major  annual  costs  of  MODEL  A  and  MODEL  B  include:   the  annual  cost 


.ng  the  auxiliary  equipment,  an  electric 


irly  cost  of  n 


of  owning  the  collector,  storagf 
and  the  like;  the  annual  cost  of  own: 
furnace,  an  electric  hot  water  he; 
pump;  the  yearly  cost  of  operatinj 

over  a  specified  number  of  years  related  to  its  life  time.  The  sum  of  these 
is  usually  taken  as  a  fixed  percentage  of  investment  each  year;  for  example, 
for  a  20  year  amortization  and  8%  interest  rate,  the  annual  cost  is  0.10185 


Operating  costs  are  primarily  due  to  power  requirements  for  the 
furnace,  the  electric  hot  water  heater,  the  air  conditioner  or  the  he 
and  for  pumping  water  and  moving  air  in  the  system,  summed  over  the  y 


of  glass  in  collectors,  or  any  other  costs  of  keeping  the  system  in  operating 

The  annual  costs  of  MODEL  A  and  MODEL  B,  in  dollars  per  year,  can  be 

For  MODEL  A 

CMA  =    (CcAc+Cst+Ce)(In'    +    (Cf+Cw+Cai) (In)    +   (PA  h+PA  c) (Cr)    +  Cm  (60) 

CMB   =    (CcAc+Cst+Ce)(In)    +    (Cf+Cw+Chp)  (I^   +    ^B.h^B.c'^r'    +  Cm      (61) 
For    the   conventional   system,    the  annual   costs  are    the   sum  of    the  following: 
the  annual   cost   of   owning   an   electric    furnace,    an   electric   hot  water  heater, 
an  air   conditioner,    and    the  annual   power   requirements   of   operating   them.      It 
is    formulated  as: 

CCON  "    <Cf«V«.i><i;>    +    (PCON,h+PCON,c><Cr>  (62) 

Since    it   has  been  assumed    that   PA>C   =  PB>C   =   P^.    the  annual   dollar 
savings   for   MODEL  A  and  MODEL   B   can  be   calculated  by: 
For   MODEL  A 

(DS)A  =   CCON-   CMA  =    (PCON,h-   PA,h)  (<V  "(CcAc+Cs t+Ce)  ( V^m  (63) 

For   MODEL  B 

(DS)B   "   CCON-    SfB   =    ^COH.h-'B.h^V-^c^Bt^tV 

"  ^p-^i^1^"^  (64) 

The  annual  dollar  savings  for  MODEL  A  and  MODEL  B  for  collector  areas  of 
400,  600,  800,  and  1000  ft2  at  electrical  rates  of  $0.02,  0.03,  0.04,  and 
0.05  per  KWh  were  calculated  and  are  presented  in  Table  6.1.   The  calculations 


L  Dollar  Savings  for  MODEL  A  and  MODEL  B  fc 
:tor  Areas  of  400,  600,  800,  and  1000  ft2  a 
Rates  of  $0.02,  0.03,  0.04,  and  0.05/KWh. 


Annual  Dollar  Savings   ($) 


The  capital  cost  per  unit  area  of  collector  ,  C  =  $8.1/ft 

The  capital  cost  of  storage'',  Cgt  =  ($0.082/lb)  (lb  water  stored) 

The  capital  cost  of  additional  equipment",  C  =  $500 

Annual  factor  of  investment,  In  and  r  =  0.10185**" 

The  cost  difference  between  a  heat  pump  and  an  air  conditioner, 

C,   -  C  .  =  $1000** 
hp    ai 

The  annual  cost  of  maintenance,  C  «  $100 

The  value  of  PA>h  and  PB>h  in  the  above  calculations  come  from  the 
appropriate  terms  of  the  total  auxiliary  energy  required  in  Table  5.3  by 

auxiliary  energy  comes  from  electricity.   The  value  of  P^  is  the  sum  c 
the  total  heating  loads  plus  6.81  million  Btu  of  hot  water  load  in  the 
summer.      The  term  was  added  since  the  solar  systems  were  assumed  to 

The  annual  dollar  savings  of  MODEL  A  and  MODEL  B  was  plotted  against 
collector  area  in  Figure  6.1.  It  can  be  seen  in  Figure  6.1  that  MODEL  B  i 
better  than  MODEL  A  for  the  costs  assumed  in  this  study  and  that  the  elect 
rate  must  be  at  about  $0,036  and  at  $0,025  per  KWh  respectively  before  any 
dollar  savings  for  MODEL  A  and  for  MODEL  B  can  occur.  The  optimum  collect 
size  for  MODEL  A  at  $0.036/KWh  is  about  600  ft2,  and  that  for  MODEL  B  at 
$0.025/KWh  is  about  500  ft2.   Since  the  current  winter  electrical  rate  for 


From  reference  (11). 
'"From  reference  (25). 

Corresponding  to  20  year  amoratization  with  8%  interest  rate. 
''(3784  Btu/hr)(12  hr/day)(30  day/month)  (5  month)  =  (6.81  x  106  Btu) 


Parameter  =  Cost  of 
Electricity,  $/KWh 


Figure  6.1  also  shows  that  the  dollar  savings  for  solar  energy  systems 

optimum.   Thus,  small  errors  in  the  estimate  of  the  long-term  thermal  per- 
formance of  solar  energy  systems  will  have  little  effect  on  the  resultant 

effect  the  results,  a  calculation  was  made  using  different  values  of  the 
annual  factor  of  investment  (In  and  r)  for  a  power  cost  of  $0.04/KWh.  The 
calculation  was  again  based  on  Eqs.  (63)  and  (64),  the  results  are  shown  in 
Figure  6.2.  In  the  figure,  the  parameter  value  of  0.149  corresponds  to  10 
year  amortization  with  8%  interest  rate  while  the  value  of  0.250  corresponds 
to  5  year  amortization  with  8%  interest  rate.  From  the  figure,  we  can  see 
that  the  assumed  annual  factor  of  investment  has  a  big  effect  upon  both  the 
amount  of  dollar  savings  and  the  optimum  collector  size.   The  figure  shows 

decreases  and  the  optimum  collector  size  also  decreases.   MODEL  B  is  still 
better  than  A. 

Although  no  calculation  was  made,  different  assumed  costs  for  the  col- 
lector and  the  storage  tank  would  also  have  a  great  effect  on  the  dollar 
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Figure   6.2     Annual  Dollar  Savings    for  MODEL  A   and   MODEL  B   for  Power 
Cost  at   $0.04/KWh. 


CHAPTER   7.      CONCLUSION 
The  economic    feasibility  of   direct   use   of    solar   energy   for    residential 

Three   solar  heating  and   cooling   systems  were  modeled  and    simulated   using 
the   1976-1977   hourly  weather  bureau   data   of   Manhattan.      The   simulation 

The   simulation  results   of   heating    (see   Tables    5.3  and   5.4)    show   that 
the  solar-assisted   heat    pump   system    (MODEL  B)    has   a  better    thermal   perfor- 
mance  than   the   conventional   solar    energy   heating   system   (MODEL  A).      The 
simulation   results   of   cooling    (see  Figure   5.29)    show    that    the  heat  pump 
operation   in  MODEL  c"    is   less   efficient    than   in  MODEL  B   ".      The  Continuous 
System  Modeling   Program   (CSMP)    which  was   used    in    this   study  has   proved    itself 
to  be  a   powerful   computer   language   in   simulating    the   solar    energy   systems. 
Although   it    uses   more   computer    time   for   compiling,    programming  using  CSMP 

Under  the  assumed  costs  for  solar  equipment  and  other  fees,  the  cost 
analysis  made  in  the  previous  chapter  (see  Table  6.1  and  Figure  6.1)  has 
resulted    the   following    conclusions: 

(1)      The   solar-assisted   heat  pump   system   is   superior    to    the   conventional 

the  amount   of  dollar   savings    for  MODEL  B    is   greater    than   that   for 
MODEL  A  at  any    electrical   rate. 
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(2)  Compared  with  the  conventional  heating  and  cooling  system  which 
uses  electricity  as  its  energy  source,  MODEL  B  with  an  assumed 
equipment  life  of  20  years  will  save  money  once  the  electrical 
rate  is  above  $0.025/KWh.  At  this  electrical  rate,  the  optimum 
collector   size  for   MODEL   B  would  be  about    500   ft2,    and    for    this 

would  be  about    51%. 

(3)  For  MODEL  A  with  an  equipment   life  of    20  years    to   compare   favorably 
with    the   conventional   system,    the   electrical   rate  must  be  about 
$0.036/KWh  before  a  dollar   savings   can  occur.      At    this    electrical 
rate,    the   optimum  collector   size  would  be  about   600   ft2,    and   for 
this   collector   size,    the   total  heating   load   supplied  by   solar 
energy  would   be  about    34%. 

(4)  The   greater   the   electrical   rate   is,    the   larger  will   be   the  optimum 
collector   size.      This   can  be   seen   from  Figure   6.1.      For    example, 
at    electrical   rate  of   $0.03/KWh,    the   figure   shows    the  optimum 
collector   size  for  MODEL  B   is   600   ft2,   but   at   $0.04/KWh,    it    shows 
the  optimum  collector   size  becomes   800   ft2. 

(5)  With   an  assumed   equipment    life   of   20  years,    MODEL  A   shows    that    it 
has   negative  amount   of   dollar   savings  while  MODEL  B   shows    that 

it  has  positive  amount  of  dollar  savings  at  the  present  winter 
electrical    rate  of    $0.030/KWh    for  Manhattan.      The   results   show 

is  400  ft2  for  MODEL  A  and  600  ft2  for  MODEL  B,  and  the  annual 
dollar  savings  is  about  $-109  for  MODEL  A  and  $156  for  MODEL  B 
(see  Figure   6.1) . 


(6)   The  assumed  equipment  life  has  a  big  effect  on  the  values  of 

dollar  savings  for  both  models.   As  the  assumed  life  decreases, 
the  value  of  the  annual  factor  of  investment  (I  and  I')  will 
increase,  thus  decreasing  the  amount  of  dollar  savings  for  both 
models.   The  calculated  results  still  show  that  MODEL  B  is  better 
than  MODEL  A  for  every  assumed  equipment  life  value  (see  Figure  6.2). 
Under  the  assumed  costs  for  solar  equipment  and  other  fees,  the  solar- 
assisted  heat  pump  system  shows  that  it  is  now  economically  feasible  compared 
to  the  conventional  system  which  uses  electricity  as  its  energy  source. 

than  a  conventional  heating  and  cooling  system  at  present,  since  most  of 
the  residences  in  the  United  States  use  much  cheaper  energy  sources,  -  i.e., 
gas  and  oil,  for  heating.   The  prospects  for  use  of  a  solar  energy  heating 
and  cooling  system  for  a  residence  are  thus  heightened  by  the  impending 
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ABSTRACT 

An  economic  feasibility  study  of  direct  use  of  solar  energy  for 
residential  heating  and  cooling  in  Manhattan,  Kansas  has  been  made  in 
this  work.   Three  solar  heating  and  cooling  systems  were  modeled  and 
simulated  on  a  digital  computer  using  the  Continuous  System  Modeling 
Program.   MODEL  A,  called  SOLAR  ENERGY  HEATING  AND  CONVENTIONAL  COOLING 
SYSTEM,  was  composed  of  a  solar  collector,  an  energy  storage  tank,  an 
electric  hot  water  heater,  an  electric  furnace  to  supplement  solar  energy 
in  the  winter,  a  conventional  air  conditioner  for  air-conditioning  in  the 
summer,  and  the  necessary  pumps,  piping,  and  controls.   MODEL  B,  called 
SOLAR-ASSISTED  HEAT  PUMP  FOR  HEATING  AND  COOLING,  was  similar  to  MODEL  A, 
with  the  exception  that  a  heat  pump  was  used  in  place  of  the  air  condi- 

an  air  conditioner  in  the  summer,  but  also  was  used  for  house  heating  in 
the  winter.   MODEL  C,  called  SOLAR-ASSISTED  HEAT  PUMP  FOR  HEATING  AND  HEAT 
PUMP  CHILLED  WATER  FOR  COOLING,  was  equiped  the  same  as  MODEL  B,  except 
the  heat  pump  cooling  operation  was  designed  to  cool  water  in  the  storage 
tank  during  the  summer  night  and  the  chilled  water  used  to  cool  the  house. 

The  simulations  were  made  for  a  typical  two  story  house  located  in 
Manhattan,  Kansas.   The  hour-by-hour  performance  of  each  model  was  studied 
using  the  1976-1977  hourly  weather  bureau  data  for  the  Manhattan  area. 
Several  collector  areas  were  tested  in  these  models  and  a  cost  analysis 

showed  that  the  heat  pump  cooling  operation  in  MODEL  C  was  less  efficient 
than  in  MODEL  B,  and  that,  for  solar  heating,  the  solar-assisted  heat  pump 
system  (MODEL  B)  had  a  better  thermal  performance  than  the  conventional 


solar  energy  heating  system  (MODEL  A).  Under  the  assumed  costs  for  solar 
equipment  and  other  fees,  the  results  of  the  cost  analysis  indicated  that 
MODEL  B  was  better  than  MODEL  A,  and  that,  MODEL  B  was  economically  compe- 
titive with  a  conventional  all  electric  system  at  the  present  power  rates. 
For  a  winter  electrical  rate  of  $0.030/KWh  in  Manhattan,  the  results  showed 
that  the  optimum  collector  size  was  400  ft  for  MODEL  A  and  600  ft  for 
MODEL  B,  and  the  annual  dollar  savings  were  about  $156  for  MODEL  B  and  an 
annual  dollar  loss  of  $109  for  MODEL  A. 


